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ABSTRACT 

The ability to measure metallicities and a-element abundances in individual red giant branch (RGB) stars 
using medium-resolution spectra (R x 6000) is a valuable tool for deciphering the nature of Milky Way dwarf 
satellites and the history of the Galactic halo. Extending such studies to more distant systems like Andromeda 
is beyond the ability of the current generation of telescopes, but by co-adding the spectra of similar stars, we 
can attain the necessary signal-to-noise ratio to make detailed abundance measurements. In this paper, we 
present a method to determine metallicities and a-element abundances using the co-addition of medium reso- 
lution spectra. We test the method of spectral co-addition using high-S/N spectra of more than 1300 RGB stars 
from Milky Way globular clusters and dwarf spheroidal galaxies obtained with the Keck II telescope/DEIMOS 
spectrographF] We group similar stars using photometric criteria and compare the weighted ensemble average 
abundances ([Fe/H], [Mg/Fe], [Si/Fe], [Ca/Fe] and [Ti/Fe]) of individual stars in each group with the mea- 
surements made on the corresponding co-added spectrum. We find a high level of agreement between the two 
methods, which permits us to apply this co-added spectra technique to more distant RGB stars, like stars in the 
M31 satellite galaxies. This paper outlines our spectral co-addition and abundance measurement methodology 
and describes the potential biases in making these measurements. 

Subject headings: galaxies: abundances — galaxies: dwarf — galaxies: evolution — galaxies: Local Group — 
galaxies: stellar content 



1. INTRODUCTION 

The Local Group is dominated by the Milky Way and An- 
dromeda galaxies, with large families of dwarf satellite galax- 
ies. The accessibility of these fantastic targets, particularly the 
dwarf spheroidals (dSphs) around the Milky Way, enables us 
to trace the dynamical and chemical properties of individual 
stars in order to investigate their formation. For decades, evi- 
dence has been growing that much of the material in the Uni- 
verse condensed into small dark matter halos at an early stage, 
and over a Hubble time, many of these halos contribute to the 
growth of massive galaxies in a "chaotic accretion" (Searle 
& Zinn 1978; White & Rees 1978; Cole et al. 2000; Die- 
mand et al. 2007). From observations, Helmi et al. (2006) 
showed that there were metal-poor stars in the Milky Way 
halo with [Fe/H] < -3.0 dex, which did not seem to exist in 
dwarf spheroidal galaxies, and therefore ruled out present day 
dSphs as the building blocks of the Milky Way halo. Recently, 
however, Kirby et al. (2008a, hereafter KGS08, 2009, 2010) 
determined the iron abundance distribution of vast majority of 
red giant branch stars (RGBs) in globular clusters (GCs) and 
ultra-faint dSphs, and they found a significant metal-poor tail 
of stars does exist in the Milky Way satellites, supporting the 
hierarchical formation of the stellar halo. 

To investigate the formation of big spiral galaxies, espe- 
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dally their stellar halo accretion history, chemical abundance 
patterns obtained from individual stars are crucial indicators 
(Wheeler et al. 1989; Worthey 1994; Mannucci et al. 2010). 
Stars produce a diversity of elements through nucleosynthe- 
sis, which are dispersed into the interstellar medium, and 
are then mixed with material in subsequent star formation. 
Generally, iron-peak elements, such as vanadium and iron, 
are mainly generated by Type la supernovae (SNe la — Tins- 
ley 1980) which are important contributors to the total iron 
fraction in galaxies (Greggio & Renzini 1983). Correspond- 
ingly, the a-elements, like oxygen, magnesium, silicon, cal- 
cium, titanium, etc., are produced in core collapse supernovae 
(SNe II) whose progenitors are massive stars with typical stel- 
lar masses greater than 9M G (Wheeler et al. 1989; Woosley 
& Weaver 1995; Gilmore 2004). Furthermore, compared to 
Type II supernovae, which have a timescale of around ~ 10 
million years (Pagel 1997; Woosley & Janka 2005), Type la 
supernovae have a longer timescale of at least ~ 1 Gyr (Mat- 
teucci & Recchi 2001; Ishigaki et al. 2012). Hence, a plot 
of [a/Fe] versus [Fe/H] tells us about the relative contribu- 
tion of SNe la and SNe II in star formation and evolution of 
a stellar system as a function of time, and can be used as a 
clock to measure the intensity of star formation at early stages. 
Additionally, comparing the distribution of metallicity and a- 
element abundances between the different stellar components 
of the Milky Way and its companions is a way of determining 
their evolutionary relationships. 

Photometry is commonly used to determine the metallici- 
ties of old RGB stars. This technique uses the locations of 
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RGB stars in the color magnitude diagram and compares them 
to empirical relations (Armandroff et al. 1993), fitting func- 
tions (Saviane et al. 2000), or theoretical stellar tracks (e.g., 
Harris & Harris 2000; Mouhcine et al. 2005; Lianou et al. 
2010) to derive photometric estimates of the metallicity. The 
metallicities of distant stars in the Andromeda system have 
also been estimated in this way. Da Costa et al. (1996, 2000, 
2002) used HST/WFPC2 to measure photometric metallici- 
ties and obtained age estimates using the RGB and horizontal 
branch in three dwarf spheroidal galaxies, Andromeda I, An- 
dromeda II, and Andromeda III, around M31. Kalirai et al. 
(2010) studied hundreds of RGB stars in six M31 dSphs, pro- 
ducing a luminosity-metallicity relation for dwarfs and com- 
pared it with that of their Milky Way counterparts. However, 
photometric metallicity estimates are often not as accurate as 
spectroscopic metallicity for two main reasons. First, the a 
abundance of star often has to be assumed. Second, the de- 
generacy of age and metallicity (Worthey 1994), even with 
resolved color-magnitude diagrams, is difficult to untangle. 
Lianou et al. (2011) compared photometric metallicity and 
spectroscopic metallicity based on the near-infrared Ca triplet 
analysis of RGBs in five Galactic dSphs, and they found that 
the agreement is good for metallicities between -2.0 dex to 
-1.5 dex but that a high fraction of intermediate-age stars 
would produce unreliable values. 

Deriving chemical abundances using spectroscopy is much 
more reliable, but also expensive in terms of data collection. 
A stellar spectrum contains a wealth of information, allow- 
ing us to derive important stellar parameters like effective 
temperature, surface gravity, metallicity and other heavy el- 
ements abundances. Several empirical calibrations can be 
employed to measure metallicity. Popular methods include 
Ca II K A3933 (Preston 1961; Zinn & West 1984; Beers et 
al. 1999) and the Ca near-infrared triplet calibration (Bica & 
Alloin 1987; Armandroff & Zinn 1988; Olszewski etal. 1991; 
Rutledge et al. 1997; Foster et al. 2010) with the prerequisite 
that [Ca/Fe] must be assumed. Unfortunately, these calibra- 
tions fail when [Fe/H] < -2.2 (Kirby et al. 2008b). 

The most reliable way to measure abundances is to use 
high-resolution spectroscopy (R > 20000). This technique 
has been used to analyze the detailed chemical abundance dis- 
tribution of the Milky Way system and its satellites (Shetrone 
et al. 1998, 2001, 2003; Venn et al. 2004; Tolstoy et al. 2009; 
Letarte et al. 2010). High-resolution spectroscopic measure- 
ments of individual stars beyond the Milky Way and its satel- 
lites, however, is extremely challenging. Moreover, high reso- 
lution spectroscopy is difficult to multiplex and often requires 
observing one star at a time. For these reasons, medium- 
resolution spectroscopy (R * 6000), for which abundances 
of some elements can still be well-measured, is an optimal 
choice for large sample stars in Milky Way satellites, and even 
for the M31 system (Guhathakurta et al. 2006; Koch et al. 
2007). 

Previous studies have targeted stars in the Milky Way satel- 
lites (e.g., Lanfranchi & Matteucci 2004; Shetrone et al. 
2009; Strigari et al. 2010). The largest sample was given by 
KGS08 and Kirby et al. (2009, 2010) who observed more than 
2500 RGB stars in Milky Way GCs and dSphs at medium- 
resolution using Keck/DEIMOS (Faber et al. 2003). They 
performed a multi-element abundance analysis using pixel- 
to-pixel matching based on a grid of synthetic stellar spec- 
tra. They verified their abundance measurement technique 
by comparing their results with the abundances derived from 
high-resolution spectra. In their study, in addition to finding a 



long metal-poor tail in ultra-faint dSphs, matching that of the 
Galactic halo, they derived trends for individual a-elements 
(Mg, Si, Ca, Ti) versus metallicity, demonstrating that there 
was not much metal enrichment before the onset of SNe I 
(Kirby et al. 2011a, 2011b). 

Observations of M31's halo present a more complex accre- 
tion history (Ibata et al. 2001; Choi et al. 2002; Reitzel & 
Guhathakurta 2002; McConnachie et al. 2004; Guhathakurta 
et al. 2006, 2010; Kalirai et al. 2009, 2010; Collins et al. 
2010; Tanaka et al. 2010), indicating that a study of abun- 
dance trends in its satellites could be more interesting and 
challenging. However, the larger distance and fainter appar- 
ent magnitudes hamper any detailed investigation of individ- 
ual stellar abundances for dwarf satellites in M31, although 
there have been attempts to obtain detailed chemical abun- 
dance patterns in M3 1 globular clusters using high-resolution, 
integrated-light spectroscopy (Colucci et al. 2009). One way 
to address this problem is by co-adding many similar stars to 
produce spectra with higher signal-to-noise (S/N). This pro- 
cess, if done properly, would extend our ability to obtain el- 
emental abundance analysis to larger distances, but there is 
also the potential for introducing biases. 

The co-addition of spectra to measure the ensemble proper- 
ties of similar objects has a long history in spectral processing. 
Adelman & Leckrone (1985) used co-addition for the ultravi- 
olet and optical region of spectrum. Holberg et al. (2003) co- 
added multiple observations of individual white dwarf stars 
to enhance the signal-to-noise ratio and combined them into 
a single spectrum. Gallazzi et al. (2008) used co-added spec- 
tra of galaxies with similar velocity dispersions, absolute r- 
band magnitude and 4000A-break values for those regions 
of parameter space where individual spectra had lower S/N. 
Most recently, Schlaufman et al. (2011, 2012) compared the 
average metallicities and a-element abundances between the 
elements of cold halo substructure (ECHOS) and the kine- 
matically smooth stellar inner halo along lines of sight in the 
Sloan Extension for Galactic Understanding and Exploration 
(SEGUE) by co-adding spectra of the metal-poor main se- 
quence turnoff stars identified by the SEGUE Stellar Param- 
eter Pipeline. Using noise-degraded spectra, they found that 
the mean square error (MSE) of abundances derived from co- 
added spectra is from 0.05 dex for metal-rich (iron-rich) stars 
to 0.2 dex for [Fe/H] and [a/Fe] for the most metal-poor stars. 
In a test using the globular clusters M13 and Ml 5, the MSE 
roughly equals to 0. 1 dex for most of the metal-rich stars and 
to 0.2 dex for metal-poor stars, in both [Fe/H] and [c/Fe]. 
However, the individual abundances of cr-elements are inac- 
cessible in their measurements. In this work, we use weighted 
spectral co-addition of RGB stars that share similar intrinsic 
properties to increase the S/N so that individual a-elements 
can be measured. 

All data used in this paper are from the same data sets used 
in KGS08 and Kirby et al. (2009, 2010), which involves thou- 
sands RGB stars in Milky Way GCs and dSphs. We test 
our co-addition method by measuring metallicity ([Fe/H]) 
and four ar-element abundances ([Mg/Fe], [Si/Fe], [Ca/Fe], 
[Ti/Fe]). We used the same definition of chemical abundances 
as defined inKGS08Q 

First, we present a short summary of the observations and 

1 In their work, 12 + log[n(Fe)/n(H)] Q = 7.52 is adopted in § 4 of KGS08, 
where n is number density (as adopted by Sneden et al. 1992). The abun- 
dances of other elements adopted solar composition from Anders & Grevesse 
(1989) in Table 8 of Kirby et al. (2010) 
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data reduction in §2. In §3 we describe the method of abun- 
dance determination using a synthetic spectral grid, as well 
as our method of co-addition: star selection and grouping, 
weighted co-adding, and abundance measurement. In §4, we 
present the comparison of weighted ensemble average of in- 
dividual abundances to abundances measured from co-added 
spectra. We also discuss discrepancies and biases between the 
two results. We summarize our work in §5. 

2. OBSERVATION AND DATA REDUCTION 

All the medium-resolution spectra of RGB stars used in this 
study are from KGS08 and Kirby et al. (2009, 2010), which 
include 2947 RGB stars in 8 Milky Way dSph galaxies and 
654 RGB stars in 14 Galactic globular clusters. The observa- 
tions were performed with DEIMOS on the Keck II telescope. 
The spectrograph configuration used the OG550 filter with the 
1200 line mmT x grating at a central wavelength of ~ 7800 A 
with a slit width of 077. The spectral resolution is ~1.2A to 
~1.3A (corresponding to a resolving power 6500 < R < 7000 
at 8500 A) with a spectral range of 6300-9 100A. Exposures 
of Kr, Ne, Ar, and Xe arc lamps were used for wavelength 
calibration and exposure of a quartz lamp provided the flat 
field. The DEIMOS data reduction pipeline developed by the 
DEEP galaxy redshift surve) ^]was used to extract one dimen- 
sional spectra (Newman et aL 2012; Cooper et al. 2012). The 
pipeline traced the edges of slits in the flat field to determine 
the CCD location of each slit. A polynomial fit to the CCD 
pixel locations of arc lamp lines provided the wavelength so- 
lution. Each exposure of stellar targets was rectified and then 
sky subtracted based on a B-spline model of the night sky 
emission lines. Then, the exposures were combined with cos- 
mic ray rejection into one two-dimensional spectrum for each 
slit. Finally, the one-dimensional stellar spectrum was ex- 
tracted from a small spatial window encompassing the light 
of star in the two-dimensional spectrum. The product of the 
pipeline was a wavelength calibrated, sky-subtracted, cosmic- 
ray-cleaned, one dimensional spectrum for each target. A hot 
star template spectrum was employed to remove the terrestrial 
atmospheric absorption introduced into the stellar spectra. In 
continuum determination, a B-spline was used to fit the "con- 
tinuum regions']^] of the spectra. Each pixel was weighted by 
its inverse variance in the fit, and the fit was performed itera- 
tively such that pixels that deviated from the fit by more than 
5cr were removed from the next iteration of the fit. For further 
details, please see KGS08. Table[T]lists all the stellar systems 
used in this study. 

3. ABUNDANCE MEASUREMENTS 

Kirby et al. presented a technique for multi-element abun- 
dance measurements of medium-resolution spectra, which en- 
able them to determine individual ff-element abundances of 
RGB stars in the Milky Way globular clusters and dwarf 
satellite galaxies. In this technique a large grid of synthetic 
spectra is used, so there is no restriction imposed on the 
metallicity range, overcoming the problems encountered by 
other methods. In brief, the photometric effective temper- 
ature (T e ff) and surface gravity (logg) are determined from 
isochrone-fitting on the color-magnitude diagram using three 

2 http://astro.berkeley.edu/ cooper/deep/spec2d/ 

The analysis pipeline used to reduce the DEIMOS data was developed at UC 
Berkeley with support from NSF grant AST-0071048. 

3 They call spectral regions with synthetic flux greater than 0.96 and a 
minimum width of 0.5 A "continuum regions". 



Table 1 

RGB spectral data sets 



Name 


Number of RGB stars 




Total 


logg< 1.4 


Globular Clusters 


NGC 6205 (Ml 3) 


68 


9 


NGC 7078 (M15) 


134 


26 


NGC 1904 (M79) 


58 


17 


NGC 2419 


95 


40 


NGC 7006 


11 


9 


NGC 7492 


20 


5 


NGC 5024 (M53)* 


49 


2 


NGC 6656 (M22)* 


50 


1 


NGC 5053* 


49 


1 


NGC 288* 


30 





NGC 5904 (M5) 


51 


6 


NGC 6838 (M71)* 


34 


1 


NGC 7089 (M2) 


44 


5 


Pal 13* 


10 





dSphs 


Canes Venatici I 


174 


36 


Draco 


298 


38 


Fornax 


675 


280 


Leo I 


813 


571 


Leo II 


258 


119 


Sculptor 


376 


153 


Sextans 


141 


16 


Ursa Minor 


212 


21 



Tj 



a 'A' These globular clusters were not used in this work because there were 
not enough RGB stars for coaddition. 

different model isochrones (Kirby et al. 2009) — Yonsei-Yale 
(Demarque et al. 2004), Victoria-Regina (VandenBerg et al. 
2006), and Padova (Girardi et al. 2002) — and an empirical 
color-based T e s (Ramirez & Melendez 2005). Then, they 
adopted the Levernberg-Marquardt algorithm (mpfit, written 
by Markwardt 2009) to find the best-fitting synthetic spectrum 
to the observed spectrum in several iterative steps by mini- 
mizing the x 1 calculated from the degraded synthetic spec- 
trum and the observed spectrum. Lastly, the stellar parame- 
ters of the best-fitting synthetic spectrum is presented as the 
observed one's. With different elemental masks, the abun- 
dances measured include iron abundance ([Fe/H]) and four 
a elements (Mg, Si, Ca and Ti). The four individual a el- 
ements were measured by considering only spectral regions 
most sensitive to the corresponding element. For example, 
the Mg wavelength mask covers about 20 neutral Mg lines in 
the DEIMOS spectral range. Even though all of the a ele- 
ments vary together in the synthetic spectrum, only Mg lines 
are used to determine [Mg/Fe]. For those are interested in 



depicted 



elemental masks, KGSOcQand Kirby et al. (2009fj< 
more details on the construction of the wavelength masks. 
They performed extensive comparisons of their medium res- 
olution results with high-resolution spectroscopic elemental 
abundances from previous studies to validate their technique. 
In this study, we inherit the idea of Kirby's technique and 
make proper modification to meet our demands to determine 
chemical abundances of more distant RGB stars beyond the 

4 Spectral mask, TABLE 2 

5 Section 4.6 and Figure 5 
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Table 2 

Atmospheric Parameter Grid. 



Quantity 


Minimum 


Maximum 


Step 


Number 


TeffCK) 


3500 


8000 


100 (J cff < 5600) 


22 








200 (r eff > 5600) 


12 


log g (cms- 2 ), T eff < 6800 


0.0 


5.0 


0.5 


11 


log g (cms- 2 ), r eff > 7000 


0.5 


5.0 


0.5 


10 


[M/H] (atmosphere) 


-4.0 


0.0 


0.5 


9 


[M/H] (spectra) 


-4.0 


0.0 


0.1 


41 


[or/Fe] 


-0.8 


+ 1.2 


0.1 


21 




Figure 1. The comparison of metallicity measured from medium-resolution spectra for individual stars between Kirby et al. and our codes, developed for co- 
added spectra. The left panel presents results for 8 GCs and the right panel for 8 dSphs. All these RGB stars meet the logg < 1.4 threshold. The key difference 
between these two measurements for individual stars is in our codes we used a fixed r c ff and log g derived from photometric estimates versus allowing T^g to 
float during the spectral fitting in Kirby et al.'s codes. The Levenberg-Marquardt algorithm is adopted to find the best-fitting spectrum with several iterative steps 
then give the abundance. The error bars shown include both the random errors in the fit and the systematic errors adopted from Kirby et al.. 



Milky Way. But firstly, we aim to test our co-add spectra 
method with the medium-resolution spectra used by Kirby et 
al., and the membership of these RGB stars has been con- 
firmed by Kirby et al.. 

3.1. Synthetic Spectral Library 

We accomplished all chemical abundance analysis on the 
basis of a large grid of synthetic spectra generated by Kirby 
et al. (KGS08, 2009, 2010). Based on ATLAS9 model atmo- 
spheres (Kurucz 1993; Sbordone et al. 2004; Sbordone 2005) 
without convective overshooting (Castelli et al. 1997, 2004; 
Castelli 2005), and a line list of atomic and molecular tran- 
sition data from the Vienna Atomic Line Database (Kupka 
et al. 1999), Kirby et al. synthesized spectra using the local 
equilibrium, plane-parallel spectrum synthesis code MOOG 
(Sneden 1973), which span the same wavelength range as 
the data (6300 to 9100A) with a resolution of 0.02 A. The 
synthetic spectra fail in modeling of the Ca II triplet, Mg I 



/18807 and the absorption lines of TiO. To avoid an unex- 
pected discontinuity they recomputed new ODFs for the new 
grid with DFSYNTHE code (Castelli 2005) and employed 
the solar composition of Anders & Grevesse (1989), but for 
Fe they used Sneden et al. (1992) (Kirby et al. 2009). To 
prevent unwittingly discarding the extremely metal-poor stars 
beyond the preliminary boundary of grid, they also expanded 
the synthetic spectral grid limit to [Fe/H]= -5.0 (Kirby et 
al. 2010). The value of [a/Fe] for the stellar model atmo- 
spheres would be different for each individual a-elements be- 
cause elements have been measured only with the spectral re- 
gions where are the most sensitive to the corresponding ele- 
ment. Therefore, an additional subgrid with the extra dimen- 
sion of a-element abundance ([a7Fe] a b un d) is also generated 
for more accurate measurement of [Mg/Fe], [Si/Fe], [Ca/Fe] 
and [Ti/Fe] at fixed [a7Fe] atm . This spectral library includes 
four dimensions: effective temperature (T e ff), surface gravity 
(logg), metallicity ([Fe/H] ), a abundance ([a/Fe] atm ) of the 
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stellar atmosphere. Table^gives the limited ranges and steps 
of these five parameters. This spectral grid is available on- 
line, and readers with interests in this grid are recommended 
to refer to more details in Kirby (201 lc). 

3.2. Individual Stellar Abundances 

To determine stellar abundances, we developed an indepen- 
dent code based on Kirby et al.'s technique with some refine- 
ments for the application to co-added spectra. For testing 
purposes, most of the medium-resolution spectra used have 
S/N > 20/pixel with mean S/N around 80/pixel. 

In spectral co-addition, an essential step is to rebin the spec- 
tra in preparation for co-adding. Our approach rebins each 
science spectrum onto a common wavelength region (6300- 
9100A with step 0.25 A). The same has been done to the 
degraded synthetic spectrum which is going to be compared 
with the rebinned science one. Considering our co-added 
spectra approach aims to be applied on the RGB stars of M31 
dwarf satellite galaxies whose spectroscopic temperatures are 
not available, we fixed the effective temperature, as well as 
the surface gravity, with the value derived from photome- 
try. When we measured the effective temperature and surface 
gravity, the Yonsei-Yale isochrones fitting was carried out on 
the CMD at an assumed age of 14 Gyr and [a/Fe] = +0.3 for 
all RGB stars, however, in Kirby et al.'s work they only es- 
timated log g by photometry. After setting the initial param- 
eters, we performed the abundance determination on the re- 
binned spectra. In order to verify that our method works well 
on individual stars, we redetermined chemical abundances of 
all RGB stars with our code. Figure [T] shows the comparison 
of [Fe/H] between Kirby's and ours. The stars in the Fig- 
ure 1 are also used for later co-addition test, but some stars 
whose spectra had insufficient S/N to measure a particular el- 
ement have been removed. The selection detailed is discussed 
in next section §3.3. 

3.3. Surface Gravity Restriction 

To start, we performed chemical abundance determination 
for more than 3600 RGB stars from 14 globular clusters and 
8 dwarf spheroidal galaxies. The stellar ages of these RGB 
stars are difficult to estimate but fortunately have only a small 
impact on the measured chemical abundances (Harris et al. 
1999; Frayn & Gilmore 2002; Lianou et al. 2010), so we as- 
sumed an age of 14 Gyr for all RGB stars (Grebel & Gallagher 
2004) and set [ff/Fe] = +0.3 empirically. Then, the effec- 
tive temperature (r e ff) and surface gravity (log g) of member 
stars were estimated by fitting Yonsei-Yale isochrones on the 
color-magnitude diagram. We then proceeded to measure the 
individual abundances as described in §3.2. We found some 
element abundances of some stars were unmeasurable, and we 
expected that low S/N is a possible reason. Thus, for the pur- 
poses of this test, stars for which we cannot measure a partic- 
ular element abundance were not used in the co-addition. Ad- 
ditionally, previous observations of RGB stars in M3 1 showed 
that only stars with M\ < -2.5 are accessible for spectroscopy. 
So we further introduced a cut in log g. Given the roughly lin- 
ear relationship between Mi and r e ff, this corresponds to a 
selection in log g. Figure 2 shows the linear relationship be- 
tween photometric log g and absolute magnitude in the /-band 
(Mi) for 7 dSphs (except Fornax, for which we do not have I- 
band data). From Figure [2]Mi = -2.5 roughly corresponds 
to a cut at log g — 1 .40. Thus, we only selected stars having 
photometric log g < 1 .40 for the co-addition. The number of 
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Figure 2. We use a linear relationship between photometric surface gravity 
(logg) and /-band absolute magnitude (Mi), which is converted from appar- 
ent magnitude with the extinction-corrected distance modulus, to select stars 
whose logg < 1.40 ( (Mi) > -2.5 ). The plot contains all dSph RGB stars 
except for Fornax for which we have no /-band data. The line is the best fit 
linear relation. Stars with log g < 1 .4 are used to produce co-added spectra. 

stars left for each dSph and GC after this selection is listed in 
Table [3] Only 8 globular clusters have enough stars for the 
following test. 



Table 3 

Star groups 



Name of Galaxies 




M t> 

l " good stars 


Globular Clusters 


NGC6205(M13) 


1 


9 


NGC7078(M15) 


4 


5 


NGC1904(M79) 


3 


5 


NGC2419 


5 


5 


NGC7006 


1 


9 


NGC7492 


1 


5 


NGC5904(M5) 


1 


6 


NGC7089(M2) 


1 


5 


dSphs 


Canes Venatici I 


4 


5 


Draco 


6 


5 


Fornax 


30 


8 


Leo I 


44 


10 


Leo II 


11 


8 


Sculptor 


16 


8 


Sextans 


2 


5 


Ursa Minor 


2 


5 



a Ntii, is the number of bins for each stellar system. 

b N goo j stars is the minimum number of good stars we set for each bin. 



3.4. Grouping and co-addition 

We consider photometric effective temperature (T e e) and 
the photometric metallicity estimate ([Fe/H] p h t), respec- 
tively, to organize the remaining stars into groups for co- 
adding. The photometric metallicities are also derived by 
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Yonsei-Yale theoretical isochrones fitting with an age of 14 
Gyr (See §3.3). We used a cut at logg < 1.40 to ensure that 
all stars lie in a range of about 1 dex in logg. Moreover, the 
synthetic spectral measurements use neutral metal lines only 
which are nearly insensitive to surface gravity. Therefore, 
log g barely changes the strength of spectral features, making 
it acceptable not to include log g in the binning. 

The goal of this study is, for each grouping, to compare 
the weighted average abundances of the individual stars (the 
input) with the abundances measured on the co-added spec- 
tra (the output). For this purpose, it is important that each 
star used has a measurable abundance. Our pipeline is able to 
measure [Fe/H] in all sample stars, but for some of them, in- 
dividual a elements (e.g., Mg) may not be measurable due to 
the quality of the spectra. Including these spectra in the group- 
ing would bias the measurement on the co-added spectrum 
(but contribute nothing to the weighted average), so we are 
very careful to construct separate groupings for each element 
measured, i.e., a co-added spectrum for testing [Mg/Fe] con- 
sists only of stars that have reliable [Mg/Fe] measurements 
individually. This allows us to use the maximum number of 
available stars for testing each element. 

After ranking member stars by their (r e ff) p h ot and [Fe/H] p h ot 
separately, we make sure that each group contains at least 5 
stars for which all five elemental abundances are measurable 
individually. We set 8 as the minimum number of stars for 
Fornax, Leo II and Sculptor, and 10 for Leo I for their large 
number of stars. Table 3 lists details for each GC and dSph. 
When we co-add spectra together for one bin, we produce 
five different co-adds, that is for each elemental abundance 
of [Fe/H], [Mg/Fe], [Si/Fe], [Ca/Fe], and [Ti/Fe], we only 
add the spectra whose elemental abundance is measurable in- 
dividually, and use that co-added spectrum to determine the 
corresponding elemental abundance. The bad spectral regions 
therefore make no contribution to the co-added spectrum for 
the element of interest, which is equivalent to the elemental 
abundance derived from the weighted ensemble average. Fig- 
ure [3] shows two binning scenarios in detail for 8 GCs and 8 
dSphs, by photometric effective temperature and metallicity. 
As the left panels show, stars having expanded distribution in 
T e g are binned more evenly, especially for Fornax, Leo I, and 
Sculptor, whose [Fe/H] p h ot are more concentrated. The im- 
pacts of this difference to the measured metallicities for these 
bins are clearly shown in the later figures. 

For the observed spectrum, we use pixel masks to remove 
bad spectral regions, like telluric absorption and cosmic rays, 
before rebinning. Keck/DEIMOS has eight CCD and the 
whole spectrum spans two CCDs, so we exclude 5 pixels near 
the end of each CCD which may cause artifacts. We then re- 
bin the spectrum onto a common wavelength range (6300A 
to 9100A) and add fluxes of the normalized rebinned spec- 
tra together weighted by the rebinned inverse variance. The 
co-addition equation is: 

_ l*j=\( x pi*el,ij/o- 2 pixel jj) 
X pixel,/ - — (1) 
^j=l 1 pixeljj 

where x pixe ijj represents the flux in the z'-th pixel of y'-th 
spectrum in a group of stars, cr pixel r is the variance of x pixe ijj, 
n represents the total number of spectra in the group, ~x p i xe i,i is 
the weighted average flux of z'-th pixel of n spectra. 

For consistency, we also create a grid of co-added synthetic 
spectra. First, for one group of stars, we pick the same number 
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Figure 3. Two scenarios have been used to bin stars. The open circles are 
individual stars used in this test. The vertical lines split the test stars into 
groups by two ways. Left panels show stars binned by photometric effec- 
tive temperature. Right panels show stars binned by photometric metallicity. 
There are five GCs (M2, M5, M13, NGC7006, and NGC7492) have one bin 
for the limited number of stars. From this comparison, the dSphs, which are 
most affected by the two scenarios, are Fornax, Leo I, and Sculptor. We detail 
the impact of different binning schemes in the discussion section. 

of synthetic spectra with same chemical abundances but with 
different and logg. The synthetic spectra are chosen to 
have the Teff and log g as determined by the photometric esti- 
mates of the observed stars. Second, we smoothed all the syn- 
thetic spectra with a Gaussian filter to match the spectral res- 
olution science spectra. Then, we co-add the synthetic spectra 
in the same way as the science spectra (Equation[T]l, with each 
pixel of the synthetic spectrum having the same weight as the 
corresponding pixel in the science spectrum. 

Figure |4]presents an example of a co-added observed spec- 
trum (blue in the top panel and black in the bottom panel and 
its best fit co-added synthetic spectrum (red in the bottom 
panel). Stars with similar stellar properties should have the 
most similar spectra which are then used to determine chem- 
ical abundances and stellar parameters. When we carried out 
the co-addition, we also avoided the calcium triplet, which 
failed in spectra synthesis. When the comparison was exe- 
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Figure 4. Top: Spectra of 5 RGB stars belonging to a single bin. These stars were selected to have logg < 1.4 and are ranked by photometric [Fe/H]. Each star's 
individual spectrum is plotted in black, and the individual abundances of stars are shown at the bottom right. The co-added spectrum for these 5 stars is plotted 
below in blue over the 5th star for comparison. When we did the co-addition process, we masked the calcium triplet and bad regions that may bias the abundance 
measurements. Bottom: the co-added spectrum (black) and its best fit co-added synthetic spectrum (red). We fixed every star's effective temperature and surface 
gravity (log g) by its corresponding photometric values, then created a grid of co-added synthetic spectra with a range in abundance values. The red spectrum 
is the best fit co-added synthetic spectrum for these 5 stars, the abundance of the best fitting synthetic spectrum and weighted average abundance are shown for 
comparison. 



cuted, only the spectral regions determined to be most sensi- 
tive to the element were considered in the abundance determi- 
nation, as showed in the bottom panel of Figure |4] 

3.5. [Fe/H] and [a/Fe] Determination 

Following KGS08 and Kirby et al. (2009), we determine 
individual chemical abundances with elemental "masks" that 
only cover the wavelength regions that are sensitive to that 
particular element. Each element, therefore, has its own mask 
but, for [a7Fe] atm the mask is a combination of the four a- 
element masks. We adopt the Levenberg-Marquardt algo- 
rithm implemented in IDL programs written by Markwardt 
(2009) to find the best match to the co-added observed spec- 
tra. This algorithm finds the best co-added synthetic spectrum 
by minimizing x 1 of pixels between the co-added observed 
spectrum and the co-added synthetic spectrum over iterations. 
The detailed procedure of our abundance determination is de- 
scribed in 1 1 steps: 

1. [Fe/H], first pass: For a bin of stars, the co-added ob- 



served spectrum is compared with the co-added syn- 
thetic spectra. In this pass only the spectral regions 
that are sensitive to Fe absorption are considered. The 
[ff/Fe] atm (total a-element abundance) and [a7Fe] a b U nd 
(abundance of individual element, Mg, Si, Ca, Ti) were 
fixed at 0, while the [Fe/H] value was varied until the 
best fit synthetic spectrum was found. When individual 
synthetic spectra were chosen, the r er i and logg were 
set to be photometric values of individual science stars 
in that group and the abundances were set to be the 
[a//Fe] atm , [a/Fe] a bund and [Fe/H] values. The derived 
[Fe/H] of the best fit synthetic spectrum was fixed and 
used in later steps. 

2. [a/Fe] atrn , first pass: All the variables were fixed except 
the total a-element abundance, [a7Fe] atm . 

3. Continuum refinement: To refine the continuum of the 
observed spectrum, the co-added science spectrum was 
divided by the co-added synthetic spectrum with the pa- 
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rameters determined in step 1 and 2, resulting in a quo- 
tient spectrum without absorption lines. Then we used 
a B-spline with a breakpoint spacing of 250 pixels to 
fit the quotient. Finally, we divided the co-added ob- 
served spectrum by the spline fit. This forces a better 
continuum match between the observed and synthetic 
spectra. 

4. Step 1-3 were repeated until [Fe/H] and [o7Fe] atm 
changed by less than 0.001 dex between consecutive 
passes. 

5. Sigma clipping: In this pass, we masked pixels whose 
absolute difference from the best-fit synthetic spectrum 
with the parameters determined in step 3 exceeded 2.5 
times their standard deviations. Then we used this 
sigma-clipped spectrum to repeat step 1. All param- 
eters except [Fe/H] were fixed. The resulting [Fe/H] 
value is the second pass value. 

6. [a7Fe] atm , second pass: We repeated step 2 with the 
sigma-clipped spectrum and the [Fe/H] determined 
from step 5. 

7. [Fe/H], third pass: Step 5 was repeated with the value 
of [a7Fe] atm determined from step 6. 

8. [Mg/Fe]: Only the spectral regions sensitive to Mg ab- 
sorption lines were considered. We fixed [Fe/H] and 
[o7Fe] atm at the values attained from steps 7 and 6, and 
varied [a7Fe] a b un d to measure [Mg/Fe]. 

9. [Si/Fe]: We repeated step 8 but replaced [Mg/Fe] with 
[Si/Fe]. 

10. [Ca/Fe]: We repeated step 8 but replaced [Mg/Fe] with 
[Ca/Fe]. 

11. [Ti/Fe]: We repeated step 8 but replaced [Mg/Fe] with 
[Ti/Fe]. 

4. COMPARISON OF CO-ADDED SPECTRAL ABUNDANCES AND 
WEIGHTED AVERAGE ABUNDANCES 

The test we conduct in this study is to measure abundances 
on co-added spectra and see if the results we get are consis- 
tent with the measured abundances of the input RGB stars. 
In future applications, the individual star will be too faint to 
measure abundance, so any biases should be anticipated by 
using this nearby star sample. In this section, we compared 
the weighted-average abundances of each bin with the abun- 
dances measured from the co-added spectra in order to test the 
feasibility of the co-addition method. The spectra we used in 
this test have a relatively high S/N, and the final goal of our 
work is to apply this method to more distant and fainter RGB 
stars beyond the Milky Way system. 

4. 1 . Weighted-Average Abundances 

In this and subsequent sections, we refer to "weighted- 
average" abundances, which are the weighted ensemble av- 
erages of the measured abundances for individual stars in a 
given bin. Correspondingly, the "co-added" abundances are 
those derived from a measurement on the co-added spectra. 

We have tried different weights to combine the individual 
abundances in bins, then compared weighted-average abun- 
dances with the abundances measured from co-added spectra. 



We found that the same weights that we used in the combin- 
ing of in dividual spectrum in co-addition (described in Sec- 
tion 3.4 1 were the best weights to use if we wanted the two 
procedures to be consistent and produce the most unbiased 
results. Taking [Fe/H] as an example, we used the elemental 
mask for Fe, which covers the wavelength regions used to de- 
termine [Fe/H], resulting in an inverse variance array. The av- 
erage inverse variance across the entire spectrum, calculated 
as in Equation[2]and Equation|3] is then used as the weight for 
that star when combining the individual abundances together 
to create the weighted-average abundance for that bin. The 
weights u>j used for individual abundances in the weighted 
average are 



mpixel 



spec,] v / ; \ypixel,i 



i,x) 



ojj(X) = 



l/cr 



spec,] 



2j=l 1 l^specj 



(2) 



(3) 



In Equation pi o~^ jiKelj j is the variance of i-th pixel of j-fh 
spectrum in the Din. mpixel is total number of pixels in a spec- 
trum. M e iementai,x is the elemental mask for measurement of X, 
where X could be abundances of [Fe/H], [Mg/Fe], [Si/Fe], 
[Ca/Fe], or [Ti/Fe]. M e i e mentai,x is a binary array in which 
only the pixels that most sensitive to corresponding element 
X absorption lines are set to 1. o 2 c ■ is the weighted variance 
for the whole j-th spectrum. In Equation|3] tOj{X) denotes the 
weight of X for j-th star, and there are n stars in that group. 
Then, the weighted average abundance is 



Xbit 



2 

spec,] 



(4) 



We calculated the weighted-average abundances of the 
four a-elements in the same way as [Fe/H] for each group. 
As mentioned before, some stars' individual element abun- 
dances were unavailable and these stars were not included in 
the co-added spectra, so we ignored them when determined 
weighted-average abundances. 

4.2. Errors 

We considered two kind of errors that contribute to the scat- 
ter in the abundance distribution: fitting error and systematic 
error. The fitting error was given by the Levenberg-Marquardt 
algorithm code. The MPFIT program determined the best-fit 
synthetic spectrum by minimizing x 1 an d g ave an estimate 
of the fitting error based on the depth of the x 2 minimum in 
parameter space. There are many other sources of error like 
the inaccuracy of atmospheric parameters, imperfect spectral 
modeling and imprecise continuum placement. We consid- 
ered all other uncertainties as systematic error. For individual 
stars, we used abundance error floors derived by Kirby et al. 
(2010) as systematic erro^] Therefore, the total error for in- 
dividual stars, (Ttotaij for j-th star, is calculated as: 



0~ roml,j(X) - 



(5) 



6 Table 5 of abundance error floors for five elements in Kirby et al. (2010). 
Data used in this work source from previous work of Kirby et al., hence, we 
expect the systematic errors are same for our measurements. 
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where o-f it j(X) is the fit error for abundance X of the j-th star, 
and cr svs j(X) is systematic error. The fit error and systematic 
error for individual stars should be uncorrelated, in which case 
the total error (cr to tal jPO) is simply the fit error (07,7 j(XJ) and 
systematic error (cr sys j(X)) added in quadrature, where X is 
either [Fe/H] or [cr/Fe], and a denotes Mg, Si, Ca, or Ti. 

For the weighted average abundance Xbm,wa> we estimated 
the variance for each bin weighted by cr 2 total j(X) which are 
same weights used for weighted-average abundances (see 
Equation [3] and Equation|4]): 



cr z /cr 4 

total J' spec,] 



cr' ,) 2 

spec,]* 



(6) 



the <Tbin,wa(X), then, is the weighted error of weighted mean 
abundance Xbi n ,wa- 

For abundances derived from the co-added spectra, we tried 
to use the same method in Kirby et al. (2010) to estimate the 
systematic errors. The distribution of the difference between 
the measured co-added values and the weighted mean values 
for same bins, divided by the expected errors, should be well 
fit by a Gaussian with unit variance, as shown in Equation [JJ 



X coadd Xfyjj 



y/(o~fit,coadd(X)) 2 + (o-bw, wa (X)) 2 + (o- S y S ,coadd(X)) 2 



1 



(7) 

where, cr fit, C oadd(X) is the fitting error of the co-added results, 
(Tbin,wa(X) is weighted mean error calculated from Equation[6j 
and cr svscoadd (X) represents the systematic error for the co- 
added results. These three types of errors are supposed to be 
independent to each other. In our case, however, the fitting er- 
rors (a fit,coadd(X)) and weighted mean errors {(Tbin,wa(X)) are 
already large enough, such that it is impossible to estimate 
systematic errors for the co-added results from Equation 17] 
Figure [5] shows the distributions of the difference, divided by 
the expected errors, and not including cr svs coadd (X), for all GC 
and dSph bins. The best-fit Gaussian is narrower than the unit 
Gaussian in all distributions, indicating that either the differ- 
ences of the two quantities in the numerator are too small, 
or the errors in the denominator are too large. However, this 
does not mean we have overestimated our uncertainties. Since 
we use the exact same stars for uncertainty estimating, they 
should contribute both to the weighted average and to the co- 
added spectral measurements. It is possible that the random 
and systematic errors in the mean abundance of each bin can- 
cel out to some degree. These could include errors such as 
those related to a spread of intrinsic abundances within a bin, 
a spread caused by spectral noise, and systematic errors re- 
sulting from T e s mismatch. Some of the errors are correlated 
between the two quantities in the numerator and cancel out 
when we take the differences. On the other hand, the co- 
addition enhances S/N of the spectra, then some flavors of 
systematic errors on the measured abundances, e.g., those re- 
sulting the photometric estimate of T e s being different from 
the true T e g, may indeed average down in the case of abun- 
dance determination from the co-added spectrum. So we ex- 
pect the systematic errors for co-added abundances should be 
smaller than the systematic error floors of individual stars. If 
the errors are correlated at some level, then there must be a 
non-negligible negative covariance term in the denominator. 
Both effects could be going on in our case. Here, we make 
a simplification for the systematic errors of the co-added re- 
sults, that we used the abundance error floors of individual 
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Figure 5. Distribution of the difference between the measurements from co- 
adds and weighted mean for 132 bins from 8 GCs and 8 dSphs divided by 
the errors of difference; systematic errors for co-adds are not included. Dis- 
tributions in left panels are stars binned by (7" c ff)phot. right panels show stars 
binned by [Fe/H] p |, ot . The solid blue curves are the best-fit Gaussian for dis- 
tributions. The dashed red curved are unit Gaussian with cr = 1. The areas of 
the unit Gaussian are normalized to the number of bins. 

stars for the co-added uncertainties and calculated them as: 



O~lotcd,coadd 



(X) = ^{cr fit , coadd (X)) 2 + (cr sys (X)) 2 



(8) 



where cr fu, C oadd(X) is the fit error of co-added abundance 
from Equation UJ <r sys (X) is the systematic error from Equa- 
tion]^] (Tt tai,coadd(X) is the total error used for co-added abun- 
dances. As Figure [5] shows, our error estimates are conser- 
vative. The true errors, accounting for covariance, must be 
slightly smaller. For the RGB stars of M31 satellites, there 
will be non-negligible random errors that result from the anal- 
ysis of low S/N spectra. Random errors often have a Gaussian 
normal distribution and contribute to the total errors in the 
measurements. We will discuss the effects of random errors 
to the error budget in future work of M3 1 . 

4.3. Comparison 

We expect that the abundances measured on the co-added 
spectra should match the weighted-average abundances of the 
stars that were used to produce the co-add. By using high 
S/N medium-resolution spectra of nearby RGB stars, we can 
know both the input and output abundances and make a ro- 
bust comparison. The co-added abundances for 8 GCs and 8 



dSphs are derived from co-added spectra as described in § 3.4 
and the weighted-average abu ndan ces are calculated based on 
ind ividu al abundances from |3.2| and combined as described 
in$U] 

For 8 GCs, the test is cleaner because we can reasonably 
assume no age, metallicity, or a abundance spread. Figure [6] 
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Figure 6. The comparison between weighted-average metallicity and metal- 
licity from co-added spectra for 8 GCs in the Milky Way. Top: the value of 
the vertical axis represents co-added results and horizontal axis is weighted- 
average results. For weighted-average abundances, we measured every star's 
abundance and averaged the abundances, weighting by the elemental mask 
weights that we used in the combining of the individual spectra. We selected 
stars with log g < 1 .4 for all 8 GCs stars and binned stars by photometric 
effective temperature (left) and photometric metallicity (right) respectively. 
Bottom: Residuals for co-added and weighted-average metallicity values ver- 
sus their average. Each GC has been denoted by an unique symbol and color. 



shows the comparisons between weighted-average metallic- 
ity and co-added spectra metallicity. Stars used here are con- 
strained by \ogg (< 1.4) and binned by photometric effective 
temperature (left figures) and photometric metallicity (right 
figures). GCs are monometallic, so the spreads in the pho- 
tometric metallicity of GCs are expected to arise from mea- 
surement uncertainty. Despite the limited number of stars in 
the GCs, the two results demonstrate a high level of agree- 
ment for these bins (top figures). The differences between the 
co-added metallicities and weighted-average metallicities are 
less than 0.06 dex for both binned scenarios as shown in the 
bottom figures. In FigurelT] the comparisons of a-elements on 
a bin basis also match weh\ except [Mg/Fe], for which the co- 
added results present relatively larger values compared to the 
weighted-average results, especially for NGC2419 as shown 
in the two bottom figures of residuals. 

For 8 dSphs, Figure [8] shows metallicity comparisons be- 
tween weighted-mean and co-added spectra measurements. 
The left side plots show bins of stars grouped by photo- 
metric T e ff and the right side plots show bins grouped by 
[Fe/H] p h ot . Each point represents one bin of stars. The two 
bottom plots present the difference of the two results versus 
their straight averages. The two results show a good match 
over —1.7 < [Fe/H] < -0.7, for which the difference less than 
0.05 dex. For [Fe/H] < -1.5 dex, the abundances from co- 
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Figure 7. Same as FigurelSlexcept for [u/Fe] instead of [Fe/H]. The differ- 
ent symbols denote the different GCs as defined in FigurelSlbut with different 
colors for the different a-elements. [Mg/Fe] in green, [Si/Fe] in red, [Ca/Fe] 
in purple, and [Ti/Fe] in blue. 

added spectra show a slight bias to be more metal-rich for 
both binning approaches, and stars binned by (r e ff) p hot show 
slightly less scatter. We use a high-order polynomial fit to 
compare the bias trends of the residuals for the two binning 
scenarios. There are total 115 bins for each scenario. For 
stars binned by photometric effective temperature, around 96 
percent of bins have residuals less then 0. 1 dex, and for pho- 
tometric metallicity, around 95 percent of bins have residuals 
less then 0.1 dex. 

Figure [9] shows the comparison of four a-elements for 8 
dSphs. Even after limiting our sample to \ogg < 1.4, i.e., the 
brighter stars were chosen, some stars still have unmeasurable 
elemental abundances and are not included in the bins. The 



plots in Figure 10 show the difference of two measurements 
as in the bottom panels in Figure [8] The scatter here is larger 
than for [Fe/H], especially for [Mg/Fe], for which the co- 
added spectra consistently scatter to higher values. The same 
trend seen in the GC stars shown in Figure [7] 

There are many factors that can cause a discrepancy be- 
tween these two mean abundances, such as S/N, Teff and log g, 
all of which may be interrelated. Since brighter RGB stars 
(with high S/N) often have lower \ogg and lower T e ff, it is 
difficult to untangle which of these stellar parameters affect- 
ing the abundance measurements most. For dSphs, there is a 
spread of ages and metallicities, which removes some of the 
degeneracy between S/N and T e g and log g, but our current 
data set still cannot sort this out. Figure 1 1 and Figure [12] 
show the impact of signal-to-noise values forGCs and dSphs7 
respectively. In the left panels the stars are binned by (7 , e fr) p hot 
and the stars are binned by [Fe/H] p h t in the right panels. The 
dashed lines show the standard deviations for the differences. 
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Figure 8. The comparison between weighted-average metallicity and metal- 
licity from co-added spectra for 8 dSphs in the Milky Way. Unique symbols 
with different colors have been used to denote the 8 dSphs in the plots. Top: 
the value of the vertical axis represents co-added results and horizontal axis 
is weighted-average results. For weighted-average abundances, we measured 
every star's abundance and averaged the abundances, weighting by the ele- 
mental mask weights which are used in the combining of individual spectrum 
in co-adding spectrum. We selected stars with log g < 1 .4 for all 8 dSphs stars 
and binned stars by photometric effective temperature (left) and photometric 
metallicity (right) respectively. Bottom: residuals between [Fe/H] coa( id and 
[Fe/H] wt . aV g. vs. the average of [Fe/Hloadd and [Fe/H] wt . aV g.. The curves in 
black are derived by polynomial fitting the residuals. 

Here we calculate the S/N of the co-added spectrum weighted 
by the inverse variance: 



<S/N coadd ) = 



SU(iAw)-(i/<^) 



(9) 



pixel, i 



where S /N coa dd is the signal-to-noise ratio for co-added spec- 
trum, (Tpixeij is the error of the z'-th pixel of co-added spec- 
trum. GCs have much higher S/N for their co-added spectra 
compared to dSphs. For dSph stars binned by (T e ff) p hot, the 
discrepancies are slightly pronounced with S/N decreasing, 
and seem to be more concentrated compared with the stars 
binned by [Fe/H] p h t in Figure 12 In general, however, the 
abundances are in good agreement for both binning methods, 
within 0.1 dex for GCs and within 0.2 dex for dSphs, with 
[Mg/Fe] having the most outliers. 

Figure [13] and Figure 14 consider the impact of <r e ff) for 
the abundance discrepancies. For each bin, We combined the 
individual r e jf weighted by the total error of the star to cal- 
culate the mean (T e ff). The left panels show stars grouped 
by (TefF)phot while the right panels show stars binned by 
[Fe/H] p hot- The dashed lines represent one standard devia- 
tion. There is no obvious trend shown in Figure 13 for GCs, 
but for dSphs, there appears a slight bias in [FeTR], where 
[Fe/H] coa dd is lower at lower (T e ff). Also, the differences in 
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Figure 9. Comparison for c-element abundances between weighted-average 
abundances and co-added spectral abundances for four tr-elements, Mg, Si, 
Ca, and Ti (top to bottom). Stars used for these plots are from the 8 dSph 
galaxies. Points in the left panels are stars binned by T c ff and right hand 
panels are binned by [Fe/HJpi,,),. The symbols are same as in Figure[8| 

[Mg/Fe] are markedly higher when 4400K < <r eff > < 4600K 
for both binning approaches. 

Figure 15 and Figure 16 consider the impact of (logg) for 
the abundance discrepancies. The photometric {\ogg) for 
each bin is derived by combining individual values weighted 
with total error of star. Like (T e tf), the differences in [Mg/Fe] 
are tied to (logg), particularly for stars binned by T e g with 
larger scatter for (logg) > 0.9. The discrepancies in metallic- 
ities and a-elements for GCs seem to be not very sensitive to 
(S/N coa dd), (7df), or (logg), and we suspect the limited num- 
ber of bins (only 17 bins for GCs) is one reason. 

4.3.1. Photometric metallicity 

In order to make a robust comparison, we also compared 
metallicity derived from co-added spectra ([Fe/H] coa dd) to the 
weighted average of photometric metallicity ([Fe/H] p hot) as 
shown in Figure 17 We combined the photometric metallici- 
ties of member stars in each bin weighted by their photometric 
errors. In the top panels of Figure 17 the photometric metal- 
licities of weighted average bins present scatter for M15 and 
NGC2419, but this is still acceptable given the uncertainty in 
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Figure 10. Comparison for tr-element abundances. The difference between 
the weighted-average abundance and co-added spectral abundance for four a- 
elements, Mg, Si, Ca, and Ti (top to bottom) versus the mean abundance. This 
is the corresponding difference plot for Figure[9] As before, points in the left 
panels are stars binned by T c ff and right hand panels are binned by [Fe/H] p i,ot- 
We also calculate the standard deviation for each residual distribution draw 
as the dashed lines. The symbols are same as in Figure[8] 

photometric measurements. NGC7006 shows the best agree- 
ment for both binning methods. For dSphs, as shown in the 
bottom panels of Figure 17 we roughly get similar results as 
in Lianou et al. (201 1), in which mean the photometric metal- 
licity estimate only has limited reliability. There is a signif- 
icant scatter in the dSphs, except for Fornax which seems to 
show a better agreement and one possible reason is the choice 
of photometric filters. In our data set, the photometry for For- 
nax is in the B and R bands. Almost all of the other galaxies 
have photometry in V/I or M/T2. The B band is more metal- 
licity sensitive than the V or M bands. Lianou et al. (2011) 
only considered V/I photometry. Therefore, it's possible a 
large amount of the uncertainty in [Fe/H] p h t comes from a 
poor choice of filters. We will discuss the impact of our age 
assumption in §4.4.4. 

4.3.2. Degradation of spectra 

The aim for developing our co-add spectral technique is 
to use it for abundance measurements of medium-resolution 
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Figure 11. The difference between abundance measurements from co-added 
spectra versus the weighted average of individual stellar abundances plotted 
versus the S/N of co-added spectra of each bin for the 8 GCs. The left pan- 
els use bins based on photometric estimates of stellar effective temperature 
(TcffVhoti while the right panel used bins based on photometric estimates of 
metallicity [Fe/H] p h t- The dashed lines mark the range, ±<r, in each panel. 
The symbols are the same as in Figure [6] The abundance differences do not 
show a strong trend with S/N. 

spectra of RGB stars in M31's dwarf satellite galaxies. The 
S/N of these stars, however, are far lower (with a typical S/N 
of < 10/pixel) than that for the stars in this test. To complete 
this test, we reran our codes with degraded spectra and com- 
pared these measurements with the weighted-average results 
derived in §4. 1 . We have randomly chosen 24 RGB stars of 
NGC185 and calculated the weighted-average S/N of these 24 
spectra with Equation|9] We reference arm_addnoise written 
by Marble (2004) to degrade the high S/N spectra. We com- 
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Figure 12. The difference between abundance measurements from co-added 
spectra versus the weighted average of individual stellar abundances plotted 
versus the S/N of co-added spectra of each bin for the 8 dSphs. The left pan- 
els use bins based on photometric estimates of stellar effective temperature 
(Teff)phot> while the right panel used bins based on photometric estimates of 
metallicity [Fe/H] p h t- The dashed lines mark the range, ±<t, in each panel. 
The symbols are the same as in Figure [8] The abundance differences do not 
show a strong trend with S/N. 



pare the errors (<xmw) of individual reduced high S/N spectra 
with the rebinned errors (o"M3i) from the weighted-average 
S/N of 24 M31 spectra, to determine the pixels whose S/Nmw 
are higher than S/Nm3i to degrade. We degrade the high S/N 
spectrum by adding new noise to the spectrum as below: 
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Figure 13. Same as Figure [TTI except the abundance differences are plotted 
versus mean of stellar effective temperature for each bin. The symbols are 
the same as in Figure]?] 
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where / i<j and / new are the fluxes of high S/N spectrum be- 
fore and after added noise, respectively. «Gaussian is Gaus- 
sian noise. cr neK is the new error for the new flux. We de- 
graded all the spectra in each bin individually before we ran 
the co-addition. We remeasure all the chemical abundances 
described in §3.4 and §3.5 but with the degraded spectra. 
Figure 18 and Figure 19 present the comparison of [Fe/H] 
and [o'/Fe] for GCs, respectively. For comparison purpose, 
we also plot the high-resolution spectroscopic metallicities of 
GCs in red. M5 gives the best match, and NGC2419 shows a 
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Figure 14. Same as Figure |12| except the abundance differences are plotted 
versus mean of stellar effective temperature of each bin. The differences tend 
to be larger and biased toward positive values at higher T c ff values especially 
for [Fe/H] and [Mg/Fe]. The symbols are the same as in FigureJI] 

discrepancy in which the coadd metallcity is less than that de- 
termined from HRS by around 0.3 dex. [o-/Fe] shows larger 
scatter compared with the [a/Fe] measured from spectra with- 
out degrading, especially for [Mg/Fe]. The comparisons for 
dSphs are shown in Figure [20] and Figure 21 The bias trends 
are larger than the trend in Figure |H]in whicn the co-added de- 
graded spectra tend to give metal-poorer results for the metal- 
rich weighted-average groups but the inverse is true for the 
metal-poor weighted-average groups. The matches are good, 
however, for both binning approaches. The [a/Fe] measured 
from degraded spectra with low S/N like in M31 satellite RGB 
stars do reflect the mean value of these bins of stars. [Ti/Fe] 
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Figure 15. Same as Figure [TTI except the abundance differences are plotted 
versus mean of stellar surface gravity for each bin. The symbols are same as 
in Figure|6] 

shows a good match while the limited neutral lines for Mg 
with such low S/N make the accuracy worse. 

We acknowledge that the procedure we are using to degrade 
the MW GC and dSph spectra does not exactly mimic the sit- 
uation with the coaddition of the low S/N M3 1 spectra. The 
MW spectra were degraded after sky subtraction and contin- 
uum division, whereas these steps may be more difficult for 
the M31 spectra given their low S/N. 

4.4. Discussion 

The discrepancy between the co-added and weighted av- 
erage abundances indicates the intrinsic degeneracy in stel- 
lar parameter space, and it is hard for us to disentangle this 
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Figure 16. Same as Figure \l2\ except the abundance differences are plotted 
versus mean of stellar surface gravity for each bin. The differences tend to be 
larger and biased toward positive values at higher log g values especially for 
[Fe/H] and [Mg/Fe]. The symbols are same as in Figure[8] 

degeneracy just using our current data sets. Uncertainties in 
the measurements also can cause discrepancy. Understanding 
causes of the bias is essential for future abundance analysis 
once we apply the co-addition technique to RGB stars in the 
Andromeda system, for this discrepancy will surely be present 
in those abundance measurements. We also rate the two bin- 
ning scenarios in this part. 

4.4.1. [Fe/H] 

Both GCs and dSphs display trends in which measured 
metallicities from co-added spectra are likely to be higher 
than the weighted-average metallicities, except for NGC7492 
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Figure 17. The comparison between weighted-average photometric metal- 
licity and metallicity from co-added spectra for 8 GCs (top) and 8 dSphs 
(bottom) in the Milky Way. The value of the vertical axis represents co-added 
results and horizontal axis is weighted-average results. For weighted-average 
abundances, we measured every star's metallicity by isochrone-fitting with 
Yonsei-Yale theoretical isochrones. We assumed an age of 14 Gyr for all 
RGB stars and [o7Fe]= +0.3. Then we averaged photometric metallicity for 
each bin weighting by photometric errors. Here all stars satisfy logg < 1.4. 
Stars are binned by photometric effective temperature (left) and photometric 
metallicity (right) respectively. The symbols are same as in Figure|6]for GCs 
and Figure [8] for dSphs. One bin of Ml 5 has been ignored in the top-right 
panel whose weighted-average photometric metallicity is less than -3.0. 

and Fornax. Globular clusters are simple stellar systems with 
old age and a single metallicity. In our GCs sample, M15, 
M79, and NGC2419 have more than one bin and these bins 
show a scatter covering an expanded range in (S/N coac jd), 
(r e ff), and (log g) in two binning approaches. The (r e ff ) shows 
more impact on the scatter, which increasing r e ff result bigger 
scatter. The rest of GCs with one bin show a good match 
between coadd and weighted-average results. The S/Ns for 
GCs are too high to assess its impact in Figure [TT[ but for 



the degraded spectra of GCs, Figure 18 demonstrates that 
even for spectra having S/N as low as 10/pixel, the coaddi- 
tion technique can still produce reliable measurements, and 
the differences between coadded metallicites and HRS metal- 
licites are acceptable. For the scatter of bins of M15, M79 
and NGC2419, we surmise the main reason is measurement 
uncertainty. 

The dSphs have more complicated and extended star for- 
mation histories with spreads in their metallicity distributions. 
Like GCs, S/N does not seem to be a dominant factor in the 
measurement bias for both high S/N coadd spectra in the top 
panels of Figure 12 and degraded coadd spectra in Figure 20 
The mean effective temperature, as (top panels) Figure T4 
shows, is tightly related to the metallicity residuals, especially 
for T e ff > 4300 . log g seems to have a slight impact on the 
residuals of stars binned by 7^ as Figure 16 shows, which 
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Figure 18. The comparison between weighted-average metallicity and 
metallicity from co-added degraded spectra for 8 GCs in the Milky 
Way. The value of the vertical axis represents co-added results and 
horizontal axis is weighted-average results. For weighted-average abun- 
dances, we measured every star's abundance and averaged the abundances, 
weighting by the elemental mask weights which we used in the com- 
bining of individual spectrum. We selected stars with logg < 1.4 
for all 8 GCs stars and binned stars by photometric effective temper- 
ature (left) and photometric metallicity (right) respectively. The sym- 
bols are same as in Figure [6] The symbols in red are the metallicities 
measured by high-resolution spectroscopy (HRS). M2 (NGC7089): Har- 
ris catalog I http://physwww.physics.mcmaster.ca/~harris/mwgc.dat I. M5 
(NGC5904), M13 (NGC6205), M15 (NGC7078), M79 (NGC1904), 
NGC2419, NGC7006, and NGC7492: Pritzl et al. (2005). Figure|6]presents 
the same information based on non-degraded spectra with the observed S/N. 




Figure 19. Same as Figure [T8] except for [ff/Fe] instead of [ Fe/H ], The 
different symbols denote the different GCs as defined in Figure [18] but with 
different colors for the different a-elements. [Mg/Fe] in green, LSi/Fe] in 
red, [Ca/Fe] in purple, and [Ti/Fe] in blue. FigurelTlpresents the same infor- 
mation based on non-degraded spectra with the observed S/N. 

just demonstrates what we note in §3.4, which is that our syn- 
thetic spectral measurements use neutral metal lines which 
are insensitive to surface gravity. The detailed comparisons 
of medium-resolution results and HRS results for individual 
stars of dSphs can be found in Kirby et al. (2010). 

Effective temperature has an important impact on the trends 
of residuals. Metallicities determined from the co-added spec- 
tra of metal-rich stars with low effective temperature (e.g., 
Fornax) tend to be more metal-poor compared with metal- 
licities derived from the weighted average method. As the 
effective temperature increases, the metallicity measurement 
residuals for co-added spectra become more metal-rich. In 
spectral synthesis, when r ef f increases, it affects [Fe/H] be- 
cause more metals are needed in the synthetic atmosphere to 
compensate for the weakening lines caused by increasing tem- 
perature (KGS08), and vice versa for metal-rich stars with 
lower effective temperature. In spectral co-addition, this ef- 
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Figure 20. The comparison between weighted-average metallicity and 
metallicity from co-added degraded spectra for 8 dSphs in the Milky Way. 
The symbols are same as in Figure [8] The value of the vertical axis repre- 
sents co-added results and horizontal axis is weighted-average results. For 
weighted-average abundances, we measured every star's abundance and av- 
eraged the abundances, weighting by the elemental mask weights which are 
used in the combining of individual spectrum in co-adding spectrum. We 
selected stars with log g < 1.4 for all 8 dSphs stars and binned stars by photo- 
metric effective temperature (left) and photometric metallicity (right) respec- 
tively. Figure|8]presents the same information based on non-degraded spectra 
with the observed S/N. 

feet seems to be amplified, probably because there is a range 
of effective temperature for each bin, and the stars with higher 
effective temperature favor the more metal-rich synthetic one 
to compensate the temperature-induced weakening lines, and 
vice versa. Schlaufman et al. (2011) noticed the effect of the 
effective temperature on the co-addition method, and they em- 
ployed 500K bins in their stellar co-additions. Our bins of 
stars do not have a fixed range in effective temperature but 
r e ff for each bin are in close proximity, a range smaller than 
500K, as shown in the left panels of Figure [3] Since T^f and 
log g are correlated, we expect log g would have the same im- 
pact on metallicity, but it shows a less prominent impact than 
Jeff on metallicity. The intrinsic metallicity dispersion of in- 
dividual systems also plays an essential role in the individ- 
ual discrepancy of different stellar systems. Schlaufman et al. 
(201 1) determined the average metallicity of M13 and M15 by 
co-addition but using metal-poor main sequence turnoff stars 
from the Milky Way halo. At lower spectral resolution, they 
measured [Fe/H]= -1.7+0.15 forM13 and [Fe/H]= -2.4+0.2 
for M15. Our results are more accurate, and the bias is less 
than 0.01 for M13 and less than 0.06 for M15 (Figure|6]l. 

4.4.2. [a/Fe] 

Compared to the iron abundance, determining individual 
tT-element abundances is more difficult, especially for dis- 
tant stars with low S/N. For both high S/N spectra and de- 
graded spectra, the comparisons of [a/Fe] between co-added 
and weighted average results show conspicuous scatter, es- 
pecially for [Mg/Fe]. For GCs, the most biased bins are 
from the metal-poor systems M15 and NGC2419, as Figure|7] 
and Figure 19 show, indicating that the measurements of a- 



elements, particularly [Mg/Fe], are very sensitive to metal- 
licity. Increasing {T e g) still has an impact on the scatter in 
multi-bin GCs. M13, with the lowest (T e g), presents an ob- 
vious bias in [Ca/Fe] and [Ti/Fe]. Compared to [Mg/Fe], 
however, these biases are still acceptable, even for the results 
derived from degraded spectra. S/N plays a more important 
role in a measurements, particular for [Mg/Fe] and [Ti/Fe] 
(Figure [TT). The same impact of S/N also shown in Figure 12 
for dSphf^ Another main reason for the scatter is the mea- 
surement uncertainty, especially for [Mg/Fe] which can be 
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Figure 21. Comparison for ff-element abundances between weighted- 
average abundances and co-added degraded spectral abundances for four a- 
elements, Mg, Si, Ca, and Ti (top to bottom). Stars used for these plots are 
from the 8 dSph galaxies. Points in the left panels are stars binned by T c ff 
and ri ght h and panels are binned by [Fe/H] p hot- The symbols are same as in 
Figure |20] Figure [9] presents the same information based on non-degraded 
spectra with the observed S/N. 

seen in the measurements of degraded spectra (Figure[19). For 
dSphs, the number of bins makes the impacts of different pa- 
rameters on the residual biases more clear. [Si/Fe], [Ca/Fe] 
and [Ti/Fe] show good agreement between co-addition and 
weighted average abundances in Figure [9] We suspect that 
the large scatter for Leo I results from the low S/N and rela- 
tively high and extended coverage in (T e g). For the outliers, 
CVnl, Sculptor, and Sextans, uncertainty in photometrically 
estimating T e g and [Fe/H], which we used to bin the stars, 
is another possible biasing source (Figure 17 1. Comparing 
Figures 12 14 and 16 it seems that S/N, T e ff, and logg are 
all possible reasons for the residual biases, especially for stars 
binned by (r e fif) p hot- Digging deeper, log g is a function of stel- 
lar mass and radius, but the mass and radius of star are also 
related to T e g. If [Fe/H] is fixed, increasing r ef f will make ab- 
sorption lines more shallow and narrow, but increasing log g 
will broaden these lines. If we consider variation in [Fe/H], 
then it becomes even more difficult to distinguish which mat- 
ters most. These small features affect the abundance measure- 



ments for both individual and co-added spectrum. Although 
stars in each bin have similar properties, the exact values of 
these important stellar parameters are different. We used the 
same chemical abundances to select the synthetic spectra for 
each bin, but the photometric and log g values used were 
those of individual observed stars. Considering the sensitivity 
of a elements to the metallicity, this assumption may play a 
role in the large scatter of [a/Fe]. Compared to other a ele- 
ments, [Mg/Fe] is the most difficult one to measure accurately 
because of weak absorption lines. Magnesium is a product of 
Type II supernova, and it is the least visible of the a elements 
in the DEIMOS spectra (KGS08). The elemental mask for 
[Mg/Fe] contains a limited number of absorption lines com- 
pared with other elemental masks, therefore results in a larger 
uncertainty for [Mg/Fe]. 

All 8 dSphs in this work have significant intrinsic spread in 
metallicity (Figure IS). The intrinsic spread in the abundance 
distribution of dwarfsatellite galaxies indicates their extended 
star formation history (Venn et al. 2004; Helmi et al. 2006; 
Cohen & Huang 2009; Kirby et al. 2009). The target selection 
of a stellar system which has radial metallicity gradient like 
Sculptor (Tolstoy et al. 2004; Walker et al. 2009; Kirby et 
al. 2009) may influence the measured abundance distribution 
for the co-added spectra method. Therefore, when we apply 
co-addition to dSphs, the effect of a metallicity spread on the 
a-elements abundance measurements with limited absorption 
lines need to be considered carefully. Particularly for the stars 
in M31 satellite galaxies with low S/N, Figures [T9| and [2T| 
show similar scatter amplitudes, demonstrating that we can 
estimate the mean values and trends for a elements from low 
S/N spectra with the coaddition technique. 

4.4.3. [a/Fe] versus [Fe/H] 

Figure [22] and Figure [23] present the individual a-element 
distribution versus metalTicity determined from co-added 
spectra for 8 GCs and 8 dSphs, respectively. The trends are 
clearly shown for the majority of bins, and only Leo I and 
Fornax display a very blurred trend in [Ca/Fe] distribution. 
We add individual stars' abundances in the figures as Fig- 
ure [24] shows. The distributions of [Ti/Fe] from Leo I and 
Fornax seem to display an upward tendency with increasing 
metallicity compared with previous results shown by individ- 
ual stars (Kirby et al. 2011b, their Figure 13). We speculate 
that it is mainly because our surface gravity restriction has re- 
moved fainter RGB stars, which are located in the metal-rich 
region in Figure [23] Yet, we still consider the co-added spec- 
tra method as an efficient and feasible tool to proceed the de- 
tailed multi-element abundance measurement with medium- 
resolution spectra. Based on the detailed or-element abun- 
dance analysis (Shetrone et al. 2001, 2003; Venn et al. 2004; 
Kirby et al. 2011b), we expect to extend our understanding 
and insight into the star formation history and galaxy evolu- 
tion to the Andromeda galaxy system and even farther sys- 
tems with this technique. 

4.4.4. Binning Scenarios 

By comparing two binning algorithms, we can assess which 
one better reflects the true values of groups of stars. The two 
parameters we used to bin, and [Fe/H] p h ot , are primarily 
estimated from isochrone-fitting on CMD. Stellar parameters 
estimated from photometry are far less accurate than the ones 
derived by spectroscopy. On the other hand, some assump- 
tions must be made to constrain the parameter space in the 



18 



Yang et al. 



[Fc/H] ^-based bins 




i : : lij.^-hascdbins 



-2.4 -2.2 -2.0 -1.8 -1.6 -1.4 -1.2 -2.4 -2.2 -2.0 -1.8 -1.6 -1.4 -1.2 
[Fe/H]„ ad [Fe/H] co „ d<1 

Figure 22. Element abundances for 8 GCs from co-added spectra. The left 
figure represents results binned by photometric effective temperature and the 
right figure is binned by photometric metallicity. 



CMD. Here we have assumed all RGB stars have an age of 
14 Gyr for both GCs and dSphs. For the dwarf spheroidal 
galaxies, this assumption is too simple because some dSphs 
are known to have extended star formation histories com- 
pared with GCs (Harris 2001; Puzia 2003; Harris et al. 2006). 
Sculptor and Sextans are dominated by old populations (Or- 
ban et al. 2008), which should show observational properties 
similar to Galactic globular clusters. However, Shetrone et al. 
(2003) proved that Sculptor has a spread in age of at least 4 
Gyr. They also demonstrated that Fornax has a complicated 
star formation history from 15 Gyr to 0.5 Gyr ago, Carina and 
Leo I have an age span of 4 - 7 Gyr and 2-7 Gyr respec- 
tively, de Boer et al. (2012) combined deep photometry with 
high resolution spectroscopy to study the detailed star forma- 
tion history and estimate ages for individual red giant branch 
stars in the Sculptor. They found that old (>10 Gyr) metal- 
poor stars, and younger, more metal-rich populations are both 
present in Sculptor. By finding a knee in the a-element dis- 
tribution at age around 10.9 Gyr they surmised that SNe la 
enrichment began (2 ± 1 Gyr) after the beginning of star for- 
mation in Sculptor. Furthermore, comparisons of the metallic- 
ities of the Milky Way dwarf spheroidal galaxies derived from 
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Figure 23. Multi-element abundances for 8 dSphs from co-added spectra. 
The left figure represents results binned by photometric effective temperature 
and the right figure is binned by photometric metallicity. The symbols are the 
same as in Figure [8] 

spectroscopy versus photometry have indicated a discrepancy 
between the two techniques (Lianou et al. 2011). Lianou et 
al. (2011) demonstrated that the assumption of old ages for 
RGB stars in dwarf satellite galaxies leads to a bias toward 
metal-poor estimates if the ages of the stars are lower, as is 
possibly the case for Fornax and Leo II. The ages of stars are 
always difficult to estimate, although one technique is based 
on using the spectroscopic metallicity, and then measuring the 
age using photometry (KGS08). 

In this work, the age affect is not enough to highly bias 
our results, despite the age-spread in dSphs. In order to test 
the impact of age, we have tried three ages in our photo- 
metric metallicity determination, 7 Gyr, 9 Gyr, and 14 Gyr, 
and we found no significant difference in these ages. Then 
we choose 14 Gyr to set to all the test RGB stars and make 
comparison of metallicity derived from co-added spectra and 
isochrone-fitting photometrically for GCs and dSphs, as Fig- 
ure \~ 



17 shows. There is a good reason, however, that [Fe/H]ph ot 
is less than [Fe/H] coa dd for Fornax. The average age is sig- 
nificantly younger than 14 Gyr, which means the RGB stars 
are bluer than the 14 Gyr isochrone. So [Fe/H] p h t will be 
too low when we assume the age is 14 Gyr. The uncer- 
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Figure 24. The same as in Figure |23| but abundances of individual stars in 
each bin have been added in (the colorful points). 

tainty in [Fe/H] p h ot may affect the results if [Fe/H] p h ot is sig- 
nificantly different from [Fe/H] determined by spectroscopy. 
However, binning by [Fe/H] p h t is still efficient, and Figure [8] 
shows that the span of metallicities derived from stars binnea 
by [Fe/H] p hot covers a larger range than the one binned by 
7^, meaning that we would most likely discover extremely 
metal-poor stars by this way in more distant stellar systems. 
These intrinsic extremely metal-poor stars play a crucial role 
in chemical evolution of these satellite galaxies and accretion 
history of galactic stellar halo (Tolstoy et al. 2001; Helmi et 
al. 2006; Kirby et al. 2008b). 

In the parameter space, the residuals of abundance binned 
by T e g span a larger range, as shown in Figure 14 and Fig- 
ure [16| illustrating that these RGB stars are very sensitive 
to !T e jf and show a larger variance with parameters such as 
S/N, r e jf and logg. However, the standard deviations of 
distributions show that chemical abundances determined by 
[Fe/H] p hot-binned stars are tighter and similar to the weighted 
average abundances. The outliers in the T e ff -binned stars are 
expected to induce the main source of the large scatter. In the 
future, we expect to combine these two parameters together, 
and to explore the impact of relationship between ages, pho- 
tometric left and logg and metallicity on binning stars. Dis- 



entangling this problem will decrease the scatter and improve 
the agreement between co-added results and weighted aver- 
age abundances. 

5. SUMMARY AND APPLICATIONS 

In this paper we present a method for spectral coaddition 
of medium resolution spectra for detailed measurements of 
multi-element abundances. Our method overcomes the low 
S/N and weak spectral features of distant RGB stars by co- 
adding spectra of similar RGB stars. We impose a surface 
gravity restriction (log g < 1 .40) for RGB stars in GCs and 
dSphs in order to enhance the quality of spectra. We group 
stars in two ways, using photometric effective temperature es- 
timates and photometric metallicity estimates, and then we 
determine detailed chemical abundances ([Fe/H], [Mg/Fe], 
[Si/Fe], [Ca/Fe], [Ti/Fe]) with elemental masks on the whole 
spectrum. We use more than 1300 well studied individual 
RGB stars from 8 GCs and 8 dSphs in the Milky Way to test 
the feasibility and accuracy of this method. This work can 
briefly be summarized as follows. 

First, we determine effective temperature, surface gravity 
and iron abundance by an isochrone-fitting technique on the 
CMD for individual RGB stars in our sample, then set these 
parameters as initial values for the following selection of syn- 
thetic spectra for each observed one and the chemical abun- 
dance comparisons. Then, individual chemical abundances 
for all RGB stars have been determined and compare well 
with previous results (KGS08; Kirby et al. 2009, 2010). 

Second, we make a cut (logg < 1.4) on candidate RGB 
stars so that our resulting sample is a good match to the M3 1 
red giant stars to which we aim to apply this spectral co- 
addition method. Then the candidate stars were binned using 
two photometric parameters, effective temperature and metal- 
licity. The science spectra belonging to one bin were added 
together, weighted by the inverse variance within specific el- 
emental masks on the pixel level. We also selected a syn- 
thetic spectrum for each observed spectrum in each bin based 
on its photometric effective temperature and photometric sur- 
face gravity. The selected synthetic spectra have the same 
metallicity and a-element abundances for each bin. Then, 
we added these synthetic candidate spectra together with the 
same weights as we used for the observed ones. 

Third, we used the Levenberg-Marquardt algorithm to se- 
lect the best fitting co-added synthetic spectrum for the co- 
added observed one by minimizing x 1 between synthetic one 
and observed one for whole spectrum on a pixel-to-pixel level 
in several iteratively steps. For each bin we determine each 
elemental abundance, [Fe/H], [Mg/Fe], [Si/Fe], [Ca/Fe], 
[Ti/Fe], separately. Only one element was considered in each 
run excluding bias probably induced by unavailable element 
abundances of individual stars. 

Fourth, we combined individual abundances in each bin 
weighted almost the same way as we performed to the co- 
added spectra, to get the weighted average abundances, and 
then we carried out comparisons with co-added abundances. 
The co-added abundances agree reasonably well with the ex- 
pected values from the weighted average abundances. 

Fifth, for (r e jf) p hot versus [Fe/H] p hot - based binning sce- 
narios, the abundance comparison and difference plots look 
similar for the two scenarios. This means the abundance er- 
rors are dominated by factors other than the spread of absorp- 
tion line strength within a bin. For metal poor stars the RGB 
isochrones are really vertical so the two binning scenarios re- 
sult in similar groups of stars. We prefer the [Fe/H] p h t based 



20 



Yang et al. 



binning scenario since our co-addition scenario for synthetic 
spectra can account for a spread in T te ff within a bin, but not 
for a spread in [Fe/H] or [ar/Fe]. 

Last, the precision of the abundances measured from co- 
added MRS does not appear to be a strong function of S/N 
over the range we have explored: co-added spectral S/N ~ 
200/pixel - 2000/pixel for GCs (see Figure [TI) and S/N 



50/pixel - 300/pixel for dSphs (see Figure nW. The S/N of 



the individual M31 RGB spectra to which we plan to apply 
this co-addition method is typically much lower than that of 
the individual spectra analyzed in this paper, but many more 
stars are co-added together in the M31 bins so that the S/N 
of the co-added M3 1 RGB spectra are expected to be compa- 
rable to that of the co-added spectra used here (Kirby et al. 
2010). The comparisons of metallicity and a elements de- 
rived from the degraded spectra with S/N comparable to RGB 
stars of M3 1 satellite galaxies and weighted average results 
demonstrate the feasibility of this technique. 

We conclude that we can safely apply this method of spec- 
tral coaddition to analyze chemical abundance patterns in 
M31 satellite galaxies using medium resolution spectra of 
RGB stars. This will be a provide a useful start for detailed 
chemical abundance exploration beyond the Milky Way sys- 
tem. 
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